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f^^ ' Abstract. We prove a theorem about elliptic operators with symmetric potential functions, defined on a function 

^\1 ^ space over a closed loop. The result is similar to a known result for a function space on an interval with Dirichlet 

, , ■ boundary conditions. These theorems provide accurate numerical methods for finding the spectra of those operators 

f*^ ' over either type of function space. As an application, we numerically compute the Morse index of constant mean 

^^ ' curvature tori of revolution in the unit 3-sphere S^, confirming that every such torus has Morse index at least five, 

^~^ ^ and showing that other known lower bounds for this Morse index are close to optimal. 

m 

rf-\ ' 1. Introduction 

^y , ', Our goal is to study the Morse index of constant mean curvature (CMC) tori of revolution in the spherical 3-space 

S^ , where the Morse index is the number of negative eigenvalues of the Jacobi operators of those surfaces. The central 
tool we use is a result about the number of nodes of eigenfunctions of those Jacobi operators. The result, proven 
J5 ' with the standard Sturm comparison technique in ordinary differential equations and closely related to classically 
known results, is proven here before being applied to the index of CMC surfaces of revolution in §^. So let us start 
by considering an operator of the form 

f^ ' on function spaces J^p over a closed loop or J-q over an interval [0, a] with Dirichlet boundary conditions: 



^ 



c° 



(N 

O^. Tp= Tp{a)= {/:M^M|/(x) = /(a: + a)}, 

2' -^0- ^o(a)- {/:[0,a]^M|/(0) = /(a) = 0}, a>0. 

^^ ■ We assume the potential function V — V{x) is real- valued and real-analytic on the closed interval [0, a], and V & Tp 

OO . when the function space J-p is used. However, we do not assume V is in J-q when the function space J-q is used, that 

^7^ ■ is, we do not assume V{Q) and V{a) are zero. 

I> . The eigenvalue problem is to find A G M and f E IFp (or Tq) that solve the second-order ordinary differential 

"/> ; equation (ODE) 

*M: (1-1) C{f)=\f. 

■ - - ' The operator —C is elliptic and it is well-known ([2], [3], [12], [19]) that the eigenvalues of C are real and form a 

discrete sequence 

Ai < A2 < A3 < ... t +00 

(each considered with multiplicity 1) whose first eigenvalue Ai is simple. The eigenvalues form a discrete spectrum, 
and corresponding eigenfunctions 

/i, /2, /s, ••■ in J^p or in J^q , C{fj) = Xjfj , 

can be chosen to form an orthonormal basis with respect to the standard L^ norm on J-'p or J-q over [0, a]. 

Let 

Ep = W/{x ^ X + a) and Sq = [0, a] 
denote the domains of the functions in J^p and J^o, respectively. The nodes of an eigenfunction f G Tp (or Tq) are 
those points of Ep (or So) at which / vanishes. When / is not identically zero, the fact that C is second-order and 
hnear implies all zeros of / are isolated and of lowest order, i.e. if /(O) = 0, then (^/)(0) 7^ 0. 

We have the following two theorems, the second of which uses a symmetry condition on V. The first theorem is 
well-known and can be proven using Sturm comparison and Courant's nodal domain theorem (see [S], [5], [7], [Sji for 
example) : 
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Theorem 1.1. Consider the operator C on the function space J-q of C°° functions over Sq with Dirichlet boundary 
conditions. Then all eigenspaces are 1- dimensional, and to find a nonzero solution f £ Tq of C{f) =^ Xf for some 
eigenvalue X, without loss of generality we may assume: 

/(0) = , (A/)(o) = i. 

Furthermore, any eigenfunction associated to the j 'th eigenvalue Xj of C has exactly j + I nodes. 

The following theorem can be similarly proven, but is a bit more complicated, because in this case the eigenvalues 
are not always simple. We will prove Theorem 11.21 here (and in the process also prove Theorem II .ip . The conclusions 
about the initial conditions in these two theorems are quite trivial; it is the conclusions about the number of nodes 
of the eigenfunctions that are of the most interest to us. 

Theorem 1.2. Consider the operator C on the function space Tp of C°° periodic functions over Sp. Suppose the 
real-analytic function V & J-p has the symmetry 

(1.2) V{x) = V{~x) \fxeR. 

Let Ai < A2 < A3 < ■•■ T +00 ^6 the spectrum of C with a corresponding basis fi, f2, fs, ... G J'p of eigenfunctions. 
Then the eigenspaces are each at most 2 -dimensional, and to find a basis for the eigenspace associated to Xj, we may 
assume: 

• When the eigenspace for Xj is 1- dimensional, we may take fj so that one of 

im - 1 , (£/.)(0) = or /,(0) = , (£/,)(0) = 1 holds. 

• When the eigenspace for Xj is 2- dimensional, and Xj = Aj+i, we may take 

im = 1 , (£/,)(0) = and /,+i(0) = , (£/,+i)(0) = 1 . 
Furthermore, any eigenfunction in J-p associated to Xj has exactly j nodes if j is even, and j — 1 nodes otherwise. 

After proving these results in Section [51 we will see in Section [3] that Theorem 11.21 gives a method to numerically 
compute the spectra of the operator C Then, in Section^ we apply that method to study the index of CMC surfaces 
of revolution in the round 3-sphere. 

In (13j . a method was given for computing the eigenvalues for the Jacobi operator of a Wente torus, involving the 
Rayleigh-Ritz method and restricting to finite dimensional subspaces of function spaces defined over tori. Then in 
[14] , both this method and a second more direct method were given for computing the first eigenvalue of the Jacobi 
operator of a Delaunay surface with respect to periodic functions, and the second method depended on Delaunay 
surfaces being surfaces of revolution. It was argued in [T3] that, although the second method was clearly the simpler 
of the two, the first method was still of value because it could compute any eigenvalue of the Jacobi operator, while 
the second method computed only the first eigenvalue. However, via Theorem 1 1.2[ the second method in [14] in fact 
extends to a method that gives any eigenvalue and hence is both simpler and equally as robust as the first method. 
Additionally, this extended second method involves only using any standard ODE solver, such as the Euler algorithm 
or the Runge-Kutta algorithm, and so has only as much numerical error as those algorithms have, whereas the first 
method involves restrictions to finite dimension subspaces for which the numerical error cannot be easily estimated 
and appears to be very much larger than for the extended second method. (This can be seen by comparing the 
respective errors of the two methods in cases where the spectra are explicitly known.) 

Certainly the first method was necessary in [13] , because Wente surfaces are not surfaces of revolution. But for 
the above reasons, the method we give here is in every way superior to the methods found in [13' and J14], in the 
case of CMC surfaces of revolution. 

A surface of revolution in the unit 3-sphere §^ is generated by revolving a given planar curve about a geodesic 
line in the geodesic plane containing this given curve. The given curve is called the profile curve and the geodesic 
line is called the axis of revolution. The profile curves of non-spherical non-flat CMC surfaces of revolution in §^ 
will periodically have minimal and maximal distances to the axis of revolution [5] . We call the points of minimal 
distance the necks, and the points of maximal distance the bulges. In general, when these surfaces close to become 
compact surfaces without boundary, they are of the following 3 types: 

• round spheres, every point of which is the same distance from a fixed point (the center) , 

• flat CMC tori, every point of which is the same distance from a closed geodesic (the axis of revolution), 

• non-flat CMC tori, where the distances from the axis of revolution to the necks and bulges are not equal. 
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Figure 1. Profile curves for surfaces that we label Ui for i — 1, 2, ..., 17 {Ui,...,Ug from left to right 
in the upper row, Uio,---,Ui7 from left to right in the lower row). These are profile curves of CMC 
tori of revolution, shown in totally geodesic hemispheres having the rotation axis as boundary. The 
images are stereographic projections from S'^ to M.^ U {oo}. The outer circle is the rotation axis, with 
profile curve inside. All of these surfaces are unduloidal, in the sense that the projections of these 
curves to the nearest points in the rotation axis are everywhere continuous local injections. 

Because these surfaces are closed, the number of negative eigenvalues of their Jacobi operators, counted with 
multiplicity and called the Morse index, is finite. The Morse index is of interest because it is a measure of the degree 
of instability of the surface. In the first two cases above, the Morse index is easily explicitly computed p^, being 1 
for the first case (this is closely related to the fact that spheres are stable jT]) and always at least 5 for the second 
case. Regarding the third case, the authors proved the following in [16 : 

Theorem 1.3. Let S be a non-flat closed CMC torus of revolution in S^ , with k bulges and k necks. Let w denote 
the wrapping number of the projection of a profile curve of S to the axis circle of revolution. Then: 

• S has index at least niax(5, 2k + I). 

• If S is nodoidal with k >2, then S has index at least niax(ll, 2fc + 5). 

• If S is unduloidal with w > 2, then S has index at least niax(6w — 1, 2A; + 4w — 3). 




Figure 2. Profile curves for surfaces that we label Ni for 1 = 1,2, ..., 11 (A^i,...,A'g from left to right 
in the upper row, TVio, iVn from left to right in the lower row). All of these surfaces are nodoidal, 
i.e. they are not unduloidal. 

The numerical results here show the lower bounds in the above theorem are very close to the true value for the 
Morse index in the case of unduloids. For example, using Table [T] and Lemnia l4.2[ the numerically computed index 
of the unduloid Ui (resp. U2, U3, ..., Un) is 6 (resp. 8, 10, 12, 14, 12, 16, 20, 24, 20, 24, 32, 28, 32, 36, 36, 44), while 
the above theorem gives the lower bound 5 (resp. 7, 9, 11, 13, 11, 15, 19, 23, 19, 23, 31, 27, 31, 35, 35, 43) for the 
index. In all cases, the lower bound in Theorem II .31 differs from the numerically computed value for the index by 
only 1, thus the lower bound is quite sharp. 

The lower bounds in Theorem 11.31 are not as sharp in the case of nodoids, but still are greater than half of the 
numerically computed value for all of the surfaces shown in Figured The numerically computed index of the nodoid 
Ni (resp. N2, N3, ..., Nil) is 12 (resp. 12, 18, 18, 24, 30, 20, 32, 34, 52, 48), while the above theorem gives the lower 
bound 11 (resp. 11, 13, 13, 15, 17, 15, 19, 21, 27, 27) for the index. 

2. Proofs of the Theorems [TTT] and [TT2I 

To prove Theorems 11.11 and 11.21 we give a series of lemmas. We first note that: 
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Figure 3. Half of each of the surfaces Ui, U2, Ur. 

• For each A, it is easily shown that the space of solutions of (jl.ip amongst functions / : M ^ M is 2-dimensional. 

• Consider the eigenvalue problem ()1.1|) over the function space J-q on the interval Sq with Dirichlct boundary 
conditions. Suppose /i and /2 are two linearly independent eigenfunctions corresponding to some eigenvalue 
A. Noting that (^/i)(0) and (^/2)(0) are both nonzero, take the linear combination f^ = (^/2)(0) • /i — 
(^/i)(a) • /2. Then /3(0) = (^/3)(0) — 0; and it follows that /a is identically zero, contradicting the linear 
independence of /i and /2. Hence the eigenvalues are simple. Hence the eigenvalues are always simple for 
the Dirichlct eigenvalue problem. Furthermore, multiplying by a scalar factor if necessary, we may assume 
the initial conditions for an eigenfunction / is /(O) — and (^/)(0) = 1. 

For the closed eigenvalue problem (|l.ip with f ^ !Fp, the eigenspace associated to any eigenvalue A is either 1 or 
2 dimensional, and we have the following lemma regarding the initial conditions to find a basis for the eigenspace: 




Figure 4. Half of each of the surfaces A^i, Nq 



Lemma 2.1. Suppose V £ J^p has the symmetry (|1.2p . Let Ai < A2 < A3 < ... t +00 ^^ ^^6 spectrum, of C 
over the space Tp with a corresponding basis /i, /2, , /a, ... of eigenfunctions. Then the eigenspaces are each at most 
2-dimensional, and to find a basis for {f € Tp\ C{f) — Xjf} for some eigenvalue Xj, we may assume: 

• When the eigenspace for Xj is \- dimensional, we may take a single eigenfunction fj e Tp such that 

either f,{0) = 1 , (£/,)(0) - or f,{0) - , (£/,)(0) = 1 . 

• When the eigenspace for Xj ~ Aj+i is 2-dimensional, we may take two eigenfunctions fj, fj+i G Tp such that 



f0) = 1 , (s/.)(o) = «"^ /.+i(o) = 



V rfx J 3 



.S^.+i)(0) = 1 • 

Proof. First we consider the case of a 1-dimensional eigenspace. Let fj G .Fp be a basis element of this eigenspace. 
If fj has neither the symmetry fj{x) = fj{—x) nor fj{—x) = —fj{x), then fj{x) and fj{—x) would be two linearly 



fj{-x) or fj{- 



-fj{x) for all 



independent eigenfunctions with eigenvalue A^, a contradiction. Hence fj{x) 

ecause multij 
Xj, we may assume either /j(0) = 1 or (^,/j")(0) = 1. Hence the first part of the lemma 



a; G M, and so fj{0) = or (^/j)(0) = 0. Furthermore, because multiplying fj by a real constant still gives a 
solution to (jl.ip with A 
is shown. 

For the case of a 2-dimensional eigenspace, any C°° solution / : M ^ M to p.ip with A 
can choose a basis fj, fj+i G J-^p with the initial conditions as in the second part of the lemma. 



Xj lies in Tp, hence we 

n 



The following lemma is known as Courant's nodal domain theorem, and the proof, which applies in our setting 
with either function space J-p or J-'q, can be found in ^ (see also ^15,). 
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Lemma 2.2. (Courant's nodal domain theorem.) The number of nodes of any eigenfunction for p.ip in J-p (resp. 
J-'o) associated to the j 'th eigenvalue Xj is at most j (resp. j + I). 



Lemma 12.21 may be strengthened using Sturm comparison, as we will see in the course of proving Theorems 11.11 
andO 

The following lemma is a slight generalization of a result in [8] : 

Lemma 2.3. Consider the following two equations 

(2.1) ^/ + (l/ + A)/ = 0, £^f + {V + \)f = 

with V as in ()1.1|) and A < A. Suppose that the first equation in (|2.ip has a solution f{x) ^ having two consecutive 
zeros at X — ^i and x — ^2, with ^1 < ^2- Let f{x) he a solution of the second equation in (j2.ip . then f{x) has at 
least one zero x = ^3 with Ci < ^3 < ^2- 

Proof. Multiplying the first equation in (|2.ip by / and the second equation in (|2.ip by /, then subtracting the first 
expression from the second and integrating, we have 

(2.2) / '(A-A)//dx- [(£/)/-(A/)/jJ' = (-|L/)(^2)/(^2)_(£/)(^,)/(^,) , 

as here f{^i) = f{£,2) = 0. Multiplying by the scalar —1 if necessary, we may assume f{x) > for x £ (^11^2), 
so (s/)(?i) > and (^/)(6) < 0. If fix) is positive everywhere in (Ci,6), then J^^{X - X)ffdx > and 

(^/)(6)/(6) - (£/)(6)/(6) < 0, contradicting (g^). Similarly, f{x) cannot be negative everywhere in (Ci,6)- 

D 

Lemma 2.4. Consider the eigenvalue problem (|l.ip on J-q over the interval Sq with Dirichlet boundary conditions, 
and with corresponding spectrum Ai < A2 < ... of simple eigenvalues. Then any eigenfunction associated with Xj has 
exactly j + I nodes. 



Proof. Denote a nonzero eigenfunction corresponding to eigenvalue Xj by fj . Lemma 12.21 implies /i has exactly two 
nodes (at a; = and x = a). Assume fj has exactly j + 1 nodal domains and let us prove /j+i has exactly j ' + 2 
nodes. From ((Tl]) we have ^/^ + Vfj + Xjfj = and J^/j+i + Vfj+i + Xj+ifj+i = 0. Let ^1,^2, ...,^j-i be the 
zeros of fj in the interval (0,a). Since Aj+i > Xj, applying Lemma |2.3[ we conclude that fj+i must vanish in each 
the intervals (0, ^1), (^1, ^2), ..., {£,j~i, a) and hence that it has at least j + 2 nodes. Lemma [2?2] implies it has exactly 
j + 2 nodes. D 

The following lemma is proven in fS]: 

Lemma 2.5. Let f and f be two linearly independent solutions of Equation (jl.ip for the same X, and suppose that 
f has two consecutive zeros ^1 and ^2 such that ^1 < ^2, then f has one and only one zero in (^1,^2). 

Proof. We may assume / is positive for all x G (Ci, ^2), then we have (^./)(^i) > and {■^f){^2) < 0. Because / and 
/ are independent, /(C,) y^O for fc- 1,2. Here £((£/)/-(|:/)/) = for all a;, so (£/)(ei)/(a) = (£/)(6)/(6). 
Hence / cannot keep a constant sign throughout the interval (Cij?2), i.e. / has at least one zero in {^i,£,2)- 

Now suppose 771 and 772 are two zeros of / in (^1, ^2). If we interchange the roles of / and / in the above argument, 
we conclude that / has at least one zero in (771,772), a contradiction. Hence / has exactly one zero in (^1,^2). D 

Lemma 2.6. Any two eigenfunctions of (jl.ip in Tp associated with equal eigenvalues have the same number of 
nodes. 

Proof. Let / and / be two eigenfunctions associated with Xj = Aj+i in the spectrum of C over J-'p. If / and / are 
linearly dependent, then the lemma clearly holds, so we assume they are linearly independent. 

Suppose / has k nodes ^1,^2, ■■■,S,k S [0,a). Then f{^e) ^ for £ = l,...,k, and by Lemma [^31 / has a unique 
node in each of (^1,^2), (^2,^3); .•.; (^fe-i,C/c) and (^fe,^i + a). Hence / has exactly k nodes. D 

Lemma 2.7. Take Xj as in Theorem (|1.2p . Let fj and fk in J-p be two eigenfunctions of C corresponding to 
eigenvalues Xj and Xk with Xj < Xk, and with either of the initial conditions as in Lemma (j2.ip . Let Uj and Uk 
denote the number of nodes in Sp of fj and fk, respectively. Lf fj and fk have the same initial conditions, resp. 
different initial conditions, then Uk > Uj + 2, resp. rik > Uj . 
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Proof. Suppose /j(0) — /fc(0) — 0, (^/j)(0) = i-sfk)iO) = 1. So fj has Uj - 1 nodes between a; = and x = a. 
Then by Lemma 12.31 fk has at least Uj nodes in the open interval (0, a), and so Uk > Uj. Since rij and n^ are both 
even, Uk > Uj + 2. 

Now suppose /j(0) — fk{0) = 1, (s/j)(0) ~ (s/fe)(0) = 0' then fj has n^ nodes fi, ...,^j in the open interval 
(0,a). Also, /j and fk have the symmetry /j(a;) = fj{—x) and /fc(a;) = fk{—x) for all a; G [0,a], by the symmetry 
(|1.2p . By Lemma [^31 /fc has a node in each interval (^^,^^+i) for -^ = 1, ...,nj — 1. Also, it has a node in {—£,l^il)^ 
so by the above symmetry, it has at least two nodes in (— Ci,'?i), implying Uk > rij and so Uk > Uj + 2. 

If fj and fk have different initial conditions, then Lemma [2T3] immediately implies Uk ^rij. D 

Lemma 2.8. Take Xj as in Theorem (|1.2p . Let fj^i, fj and fj^i inJ-p be three consecutive eigenf unctions associated 
with Xj-i, Xj and A^+i, respectively, each with either of the initial conditions as in Lemma (j2.ip . Let Uj-i, Uj and 
rij+i denote the number of nodes of fj-i, fj and fj+i respectively. Then Uj^i > Uj-i + 2. 

Proof. Here each of the eigcnfunctions fj-i, fj and fj+i has either of the two initial conditions given in Lemma 
(|2.ip . Thus two of these functions will have the same initial conditions, hence by Lemma [2?7l we have at least one of 



nj > Mj-i + 2 or 7ij+i > nj^i + 2 or rij+i > nj + 2. Lemma 12.71 also implies Uj^i < Uj < tIj+i, hence the result is 
shown. D 

All of the lemmas in this section (Section ^ immediately imply Theorems 11.11 and 11.21 

3. Computation of the spectrum of C over Tp with symmetric V 

A numerical method for computing the spectrum of C on the function space J-p is as follows: 

(1) The eigcnfunctions are in J-p and so are periodic, and the real-analytic V S J-p is assumed to have the 
symmetry ()1.2p . Theorem 11.21 implies we can numerically solve (jl.ip for / with the initial conditions just 
either /(O) = 1, (a|/)(0) = or /(O) = 0, (a|/)(0) = 1 by a numerical ODE solver, and search for the 
values of A that give periodic solutions /, i.e. give / € J-'p. Such values of A are amongst the Xj. 

(2) By Theorem 11.21 we know the eigenspaces are at most 2-dimensional. If, for some A = A^ , one of the two 
types of initial conditions in Theorem 11.21 for / gives a solution f <E J^p and the other does not, then the 
eigenspace of Xj is 1-dimensional; if both types of initial conditions give solutions / G J-'p, then the eigenspace 
of Xj is 2-dimensional. 

(3) From Theorem II. 2 1 we know that any eigenfunction / associated to Xj has exactly j nodes when j is even, 
and J — 1 nodes otherwise. So the value of j is determined simply by counting the number of nodes of /. 
Because we can determine j, we will know when we have found all Xj < M for any given M G M. 

4. AppHcation to CMC surfaces of revolution in S'^ 

As an application of the numerical approach described in Section [31 we consider CMC tori of revolution in the 
unit 3-sphere S^ C W^ and compute the spectra of their Jacobi operators. This gives us a numerical evaluation of 
the Morse index of these surfaces. 

Let S{x,y) : T — {{x,y) G M^|(x,y) = {x -\- a,y) = {x,y -\- 27r)} ^ S'^ be a conformal immersion from the torus 
T to S'^, with mean curvature H and Gauss curvature K . When H is constant, S is critical for a variation problem 
whose associated Jacobi operator is 

-A - y with V = A + 4ij2 - 2K , 

where A is the Laplace-Beltrami operator of the induced metric ds'^ = g{dx^ -\- dy^) for some smooth function 
g = g{x,y) {g is in fact real-analytic in the application here). We take iS to be a non-flat CMC torus of revolution. 
Let us define 

C = -gA-gV^£--^-V, 

where V = gV. Then the eigenvalues of C form a discrete sequence whose corresponding eigcnfunctions can be 
chosen to form an orthonormal basis for the L^ norm over T with respect to the Euclidean metric dx'^ -\- dy^. Let 

Ai < A2 < A3 < ... t -t-00 

be the spectrum of C 

By using Rayleigh quotient characterizations for eigenvalues it can be shown that C and —A — V will give the 
same number of negative eigenvalues (counted with multiplicity), although these two operators will have different 
eigenvalues. Hence we can use either £ or —A — ^ to find the Morse index of the surface S: 
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surface 


nonpositivc eigenvalues Ai o- ^2^0' ■■■> -^A;,0 (where A^ q — and A^_|_2^ q > ^)' ^°^ ^^e operator C 


t/i 


-1.28, 


-1, -1, -0.25, 


C/2 


-1.08, 


-1, -1, -0.76, -0.76, -0.51, 


Us 


-1.04, 


-1, -1, -0.87, -0.87, -0.67, -0.67, -0.52, 


Ui 


-1.03. 


-1, -1, -0.91, -0.91, -0.78, -0.78, -0.6, -0.6, -0.48, 


u^ 


-1.02, 


-1, -1, -0.94, -0.94, -0.84, -0.84, -0.714, -0.714, -0.57,-0.57,-0.48, 


Ue 


-1.64, 


-1.48, -1.48, -1, -1,-0.36, 


Uj 


-1.13, 


-1.1, -1.1, -1, -1, -0.84, -0.84, -0.64, -0.64, -0.5, 


Us 


-1.05, 


-1.04, -1.04, -1, -1, -0.94, -0.94, -0.86, -0.86, -0.77, -0.77, -0.68, -0.68, -0.64, 


U9 


-1.029 


, -1.022, -1.022, -1, -1, -0.96, -0.96, -0.92, -0.92, -0.86, -0.86, -0.79, -0.79, -0.72, -0.72, -0.66, -0.66, -0.63, 


Uio 


-1.28, 


-1.25, -1.25, -1.15, -1.15, -1, -1, -0.83, -0.83, -0.72, 


Uii 


-1.11, 


-1.1, -1.1, -1.06, -1.06, -1, -1, -0.92, -0.92, -0.83, -0.83, -0.75, -0.75, -0.71, 


Ul2 


-1.04. 


-1.03, -1.03, -1.02, -1.02, -1, -1, -.97, -.97, -0.94, -0.94, -0.896, -0.896, -0.85, -0.85, -0.8, -0.8, -0.76, -0.76, -0.73, -0.73, -0.72, 


Ul3 


-1.14. 


-1.13, -1.13, -1.1, -1.1, -1.06, -1.06, -1, -1, -0.94, -0.94, -0.88, -0.88, -0.86, 


Uii 


-1.14. 


-1.13, -1.13, -1.1, -1.1, -1.06, -1.06, -1, -1, -0.93, -0.93, -0.84, -0.84, -0.75, -0.75, -0.67, -0.67, -0.63, 


Ui, 


-1.07. 


-1.06, -1.06, -1.05, -1.05, -1.03, -1.03, -1, -1, -0.96, -0.96, -0.92, -0.92, -0.87, -0.87, -0.83, -0.83, -0.78, -0.78, -0.75, -0.75, -0.74, 


Uie 


-1.11, 


-1.1, -1.1, -1.09, -1.09, -1.07, -1.07, -1.04, -1.04, -1, -1, -0.96,-0.96, -0.93, -0.93, -0.9, -0.9, -0.89, 


Un 


-1.26, 


-1.25, -1.25, -1.23, -1.23, -1.19, -1.19, -1.14, -1.14, -1.08, -1.08, -1, -1, -0.91, -0.91, -0.81, -0.81, -0.71, -0.71, -0.62, -0.62, -0.59, 


Ni 


-1.26, 


-1.19, -1.19, -1, -1, -0.85, 


iVa 


-1.42, 


-1.31, -1.31, -1, -1, -0.696, 


Ns 


-1.43, 


-1.37, -1.37, -1.22, -1.22, -1, -1, -0.85, 


Ni 


-1.85, 


-1.76, -1.76, -1.47, -1.47, -1, -1, -0.55, 


Ns 


-1.67, 


-1.62, -1.62, -1.49, -1.49, -1.27, -1.27, -1, -1, -0.83, 


Ne 


-1.7, - 


1.67, -1.67, -1.58, -1.58, -1.43, -1.43, -1.22, -1.22, -1, -1, -0.88, 


N-r 


-1.09, 


-1.08, -1.08, -1.05, -1.05, -1, -1, -0.95, -0.95, -0.93, 


Ns 


-1.47, 


-1.45, -1.45, -1.39, -1.39, -1.29, -1.29, -1.16, -1.16, -1, -1, -0.84, -0.84, -0.75, 


Ns 


-1.18, 


-1.176, -1.176, -1.15, -1.15, -1.11, -1.11, -1.06, -1.06, -1, -1, -0.94, -0.94, -0.89, -0.89, -0.87, 


Nio 


-1.31, 


-1.3, -1.3, -1.29, -1.29, -1.26, -1.26, -1.22, -1.22, -1.18, -1.18, -1.12, -1.12, -1.06, -1.06, -1, -1, -0.94, -0.94, -0.9, -0.9, -0.89, 


Nil 


-1.19, 


-1.18, -1.18, -1.17, -1.17, -1.15, -1.15, -1.12, -1.12, -1.09, -1.09, -1.05, -1.05, -1, -1, -0.95, -0.95, -0.91, -0.91, -0.89, -0.89, -0.88, 



Table 1. Numerical estimates for the nonpositive eigenvalues of the operator C for the specific 
examples of CMC non-flat tori of revolution shown in Figures [T] and [5] The values are rounded off 
to the nearest hundredth or thousandth. 



Definition 4.1. The Morse index Ind(iS) of S is the sum of multiplicities of the negative eigenvalues of —A — V with 
function space the smooth functions from T to M. Equivalently, it is the sum of the multiplicities of the negative 
eigenvalues of C over the same function space. 
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Figure 5. Eigenfunctions associated to the eigenvalues Ai,o, •■•, A4,o, A5,o = of the surface Ui. 
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Figure 6. Eigenfunctions associated to the eigenvalues Ai,o, •■•, Ae.o, M,o = of the surface 1/2- 
A C°° function / = /(a;, y) -.T ^rM. can be decomposed into a series of spherical harmonics as follows: 
(4.1) / = X! "jm(^) cos(jy) + Uj^2{x) sin(jj/) , 

where Uj^i,Uj^2 S J-p for Sp with the given a > 0. The operator C on the function space J-'p is defined by 



C 



dx'^ 



V , 



and the spectrum 



Ai,o < A2,o < A3,o < ... T +00 
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Table 2. Here Bi is the number of eigenvalues less than —4, B2 is the number of eigenvalues in 
[—4, —1), and B^ is the number of eigenvalues in (—1, 0), all counted with multiplicity. 
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Figure 7. Eigenfunctions associated to the eigenvalues Ai^o, ••■, As.o, Ag^o = of the surface U3 



of C has all the analogous properties as those of the spectrum for C. Furthermore, by uniqueness of the spherical 
harmonics decomposition, / is an eigenfunction of C for the eigenvalue A if and only if each Uj^k, k — 1,2, is an 
eigenfunction of C for the eigenvalue A — j^. And if / is not identically zero, then some Uj^k will also be not identically 
zero. Thus we can say: 

• A is an eigenvalue for the operator C if and only if A — n^ is an eigenvalue for the operator jC for some 
neNU{0}. 
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Figure 8. Eigenfunctions associated to the eigenvalues Ai,o, •••, Aio,o, An^o = of the surface C/4. 




Figure 9. Eigenfunctions associated to the eigenvalues Ai^o, •••, ^12,0, Aia^o = of the surface U5 




Figure 10. Eigenfunctions associated to the eigenvalues Ai,o, •••, ^6,0, -^7,0 = of the surface Uq. 



if A e [-P,-{i- 1)2) for i eN 



• For any eigenvalue Xj^ < — n^, n £ N fl [2, 00), of C, with associated eigenfunction fj G J-p, the eigenvalues 
of C associated to the eigenfunctions fj ■ cos{ky) and fj ■ sm{ky), for integers k < [0, n], will be negative. 
Furthermore, we can conclude the following: 

Lemma 4.2. We have 

Ind{S) ^ Y^ Xj,o ■ £{Xj.o) , where £(A) = \ 2i - l If Xe \-P -(i - 1)^) 
jeN ^ [ T \ J J 

Proof. Let E{X), resp. Eq{X), denote the eigenspace of solutions of £(/) = A/ for smooth / : T ^ R, resp. 
£(/) = A/ for / € Tp. Then dim_E(A) == 0, resp. dimi?o(A) = whenever A is not an eigenvalue of C, resp. C, and 
is a positive integer otherwise. Then, by the uniqueness of the spherical harmonics decomposition, 

^dim^(A) = ^ I dim£;o(A)+2^dim£;o(A-j2) J ^ I ^dim£^o(A) ) +2^ ^ dimEa{X) . 

A<0 A<0 y j>l J \\<0 ) ]>lX<-p 

n 
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Figure 11. Eigenfunctions associated to the eigenvalues Ai^o, •■■, '^10,0, An^o = of the surface C/7. 



Here 5 is a CMC surface of revolution, so, following [16], we can consider 



C--K 



dip 



2v^ - 32s2t\~2 and C = 



dx' 



- 2v^ - 32sh'^v-^ 



where s G M^, t G {—s,s) \ {0}, 7 G (0, 7r/4], and s, t, 7 (we note that cot(27) is the mean curvature of S) satisfy 
the conditions (s + t)^ — Ast sin^ 7 = 1/4 and 

ste (-(16sin2 7)-i,0)u(0, (16cos2 7)-i) , 

and V — v{x) is the elliptic function 

V — 7- — - with period xq = — j — , = 

dn,(2sa:) ^ s io ^J {I ~ Q^){1 - t^ q'') 



, where r = -y/l — t^/s^ 



When st > 0, we have unduloidal surfaces. When st < 0, we have either nodoidal or unduloidal surfaces (see [IB]). 

Using the method in Section [3l we can numerically compute the negative eigenvalues of C, and can then apply 
Lemma 221 to find Ind(5). We do this for the CMC tori of revolution shown in Figures [1] and [H In [TB], it is shown 
that is an eigenvalue of C, and —1 is an eigenvalue of C with multiplicity 2. Since i{X) is discontinuous at A = and 
A = —1, it is crucial to know that both and —1 are eigenvalues of C in order to determine Ind(iS). Furthermore, 
as the eigenvalue —1 has multiplicity 2 and Ai.o must be simple, we have Ai.o < —1- (In the numerical experiments 
here, we find that is always a simple eigenvalue.) 

Tables [1] [2] and Figures [5ll32l show results of our numerical computations. 




Figure 12. Eigenfunctions associated to the eigenvalues Ai.o, •■■, A14.0, Ais.o = of the surface Us- 
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Figure 13. Eigcnfunctions associated to the eigenvalues Ai^Oj ■■■, ^^18,0, Aig^o = of the surface C/g 
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Figure 14. Eigenfunctions associated to the eigenvalues Ai^o, ■•■, Aio,o, Mi,o = of the surface C/i, 




Figure 15. Eigenfunctions associated to the eigenvalues Ai^o, •••, Ai4^o, Ai5^o = of the surface Uu. 
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Figure 16. Eigenfunctions associated to the eigenvalues A 
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Figure 17. Eigenfunctions associated to the eigenvalues Ai 
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Figure 18. Eigcnfunctions associated to the eigenvalues Ai^o, •••, Msfi, -^19,0 = of the surface U14. 




Figure 19. Eigcnfunctions associated to the eigenvalues Ai^o, •••, ^22,0, ^23,0 = of the surface U15 
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Figure 20. Eigenfunctions associated to the eigenvalues Aio, •••, Ais.o, Aig.o = of the surface U] 
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Figure 21. Eigenfunctions associated to the eigenvalues Ai.o, •.., A22.O1 A23.0 — of the surface U17. 
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Figure 22. Eigenfunctions associated to the eigenvalues Ai^o, •••, ^e.Oi A7.0 = of the surface Ni. 
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Figure 23. Eigenfunctions associated to the eigenvalues Ai,o, •••, ^6,0, ^7,0 = of the surface N2- 
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Figure 24. Eigenfunctions associated to the eigenvalues Ai.o, .-., Ago, Ago = of the surface A^3 
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Figure 25. Eigenfunctions associated to the eigenvalues Ai.o, •••, Ag.o, Ag.o = of the surface A'4. 
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Figure 26. Eigenfunctions associated to the eigenvalues Ai_o, •••, Aio,o, An^o = of the surface N5 
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Figure 27. Eigenfunctions associated to the eigenvalues Ai.o, •••, A12.0, A13.0 = of the surface Nq. 




Figure 28. Eigenfunctions associated to the eigenvalues Ai.o, •••, Aio.o, An.o = of the surface iVy. 




Figure 29. Eigenfunctions associated to the eigenvalues Ai_0) ■■■, ^4,0, Ai5^o = of the surface iVg. 
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Figure 30. Eigenfunctions associated to the eigenvalues Ai^o, •••, Aie.o, Ai7,o = of the surface Ng. 
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Figure 31. Eigenfunctions associated to the eigenvalues Ai.o, •••, A22.OJ A23.0 = of the surface N] 
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Figure 32. Eigenfunctions associated to the eigenvalues Ai,o, •••, A22.0, A23.0 = of the surface A^n. 



